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A search for the decay of the r lepton to seven charged pions and one or zero 7r° mesons was 
performed using the BABAR detector at the PEP-II asymmetric-energy e + e~ collider. The analysis 
uses 232.2 fb _1 of data at center-of-mass energies on or near the T(4S) resonance. We observe 
24 events with an expected background of 21.6 ± 1.3 events. Without evidence for a signal, we 
calculate an upper limit of B(t~ — > 47r~37r + (7r°)z/ T ) < 3.0 x 10 -7 at 90% confidence level. This is 
an improvement by nearly an order of magnitude over the previously established limit. In addition, 
we set upper limits for the exclusive decays r _ — > 47r~37r + ^ T and r _ — > 4ir~ 3tt + tt° v T . 

PACS numbers: 13.35.Dx, 13.66.De, 13.85.Rm 



The decay of the r lepton to seven charged parti- 
cles has not been observed to date. An upper limit of 
B(t~ -> 47r-37r+( 7 r ( > T ) < 2.4 x 1(T 6 at 90% confidence 
level has been set by the CLEO Collaboration The- 
oretical calculations using an effective chiral Lagrangian 
to estimate the matrix elements show that the r decay 
rate to seven charged particles is much smaller than the 
decay into five charged particles due to its smaller phase 
space. This leads to a theoretical branching fraction es- 
timate of the order of 1CP 11 , which could be enhanced 
by up to an order of magnitude if this decay proceeds via 
resonances 0. 

This analysis is based on data recorded by the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage ring operated at the Stanford Linear Acceler- 
ator Center. The data sample consists of 232.2 fb _1 
recorded at center-of-mass (CM) energies of 10.58 GeV 
and 10.54 GeV. With an expected cross section for 
t pairs for these CM energies of er TT = (0.89 ± 0.02) 
nb the number of produced r pair events N TT is 
(206.6 ±5.2) x 10 6 . 

The BABAR detector is described in detail in Ref. Q 
and only a brief description is given here. Charged- 
particle momenta are measured with a 5-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH) inside a 1.5-T solenoidal magnetic field. 
A calorimeter (EMC) consisting of 6580 CsI(Tl) crys- 
tals is used to measure electromagnetic energy. A ring- 
imaging Cherenkov detector is used to identify charged 
hadrons, in combination with ionization energy loss mea- 
surements in the SVT and the DCH. Muons are identified 
by an instrumented magnetic- flux return (IFR) . Particle 
four-momenta are reconstructed in the laboratory frame 
and then boosted to the e + e~ center-of-mass (CM) frame 
using the measured beam energies. 

We use Monte Carlo simulation techniques to estimate 
our signal efficiencies and background contamination 
from other r-pair events. The production of r pairs is 
simulated with the KK generator and non-signal r 
lepton decays are modeled with TAU0LA |(| according 
to measured rates 0. Signal events were generated 
with uniform density throughout the available phase 
space. The simulation of the BABAR detector is based 
on GEANT 4 Q. 



Events with eight charged tracks and a net charge 
of zero are selected. Tracks are required to have a 
distance of closest approach to the interaction point 
in the plane transverse to the beam axis (DOCAxy) 
of less than 1.5 cm, and a distance of closest approach 
along the beam direction of less than 10 cm. It is 
required that at least five of these eight tracks have 
a minimum transverse momentum (j>t) of 100 MeV/c, 
and 12 or more DCH hits. Photons are reconstructed 
from EMC clusters and are required to have a mini- 
mum energy of 70 MeV, more than three crystal hits, 
and a lateral energy profile consistent with that of a 
photon. The event is divided into hemispheres by a 
plane perpendicular to the thrust axis 0, where the 
thrust is calculated from all charged tracks and all 
photons in the event. The event thrust has to be larger 
than 0.9. We require all events to have one track in 
one hemisphere (the "tag-side" hemisphere) and seven 
tracks in the other hemisphere (the "signal side" hemi- 
sphere). The above requirements define the 1-7 topology. 

Monte Carlo studies have shown that the background 
from r-pair events stems from photon conversions in the 
detector material. To reduce this contribution and to ex- 
clude e + e~ — > qq (q — {u,d, s,c,b}) decays containing 
kaons we apply particle identification on the signal side. 
We demand that at least six of the tracks are identified 
as pions with high probability. We apply looser identifi- 
cation criteria to the seventh track. 

To further reduce photon conversions, we require all 
seven tracks on the signal side to have a minimum 
transverse momentum of 100 MeV/c and the ratio of 
DOCAxy /pt to be less than 0.7cm/( GeV/c). 

Further suppression of e + e~ — * qq events is achieved 
by requiring the tag side to satisfy one of the following 
criteria: 1) a tightly identified electron or muon with 
no more than one additional photon on the tag-side; 

2) a tight lepton veto and no additional photon; or 

3) a reconstructed p meson with an invariant mass of 
0.650 < nip < 0.875 GeV/c 2 , derived from the combi- 
nation of the 1-prong track with a reconstructed ir° 
candidate. The ir° candidates are formed by combining 
two photons and requiring an invariant mass between 
0.113 and 0.155 GeV/c 2 . 

After applying these requirements we calculate the 
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TABLE I: Cumulative signal efficiencies in % and r back- 
ground contributions as number of events scaled to the data 
luminosity after the various selection criteria. ID denotes 
Identification, Conv. Conversion, and Gen. Generic. 





Efficiency [%] 


Number of events 




7irv T 


77T7T J h> T 


Gen. r 


5tvv t 


5irn v T 


1-7 Topology 


23.6±1.4 


22.1T1.4 


767 


198 


187 


Particle ID 


20.7±1.3 


19.6±1.3 


108 


64 


75 


Conv. Veto 


15.8±1.0 


14.9±1.0 





4.7 


9.2 


1-Prong Tag 


10.2±0.7 


9.6±0.7 





1.7 


4.2 


Signal Region 


9.4±0.6 


9.3T0.6 





0.4 


0.8 



pseudo mass (m) of the r lepton : 

2 (Eb eam — Eit T )(E 77T — P^) 



m 



(1) 



The pseudo mass is an approximation of the invariant 
mass of the tau, where the neutrino's flight direction is 
approximated by the combined momentum vector 
of the seven charged tracks, and its energy is taken to 
be the difference between the beam energy E^ eam in the 
center-of-mass system and the combined energy Ej^ of 
the charged tracks. The pseudo mass allows for a better 
discrimination between signal events and background 
from e + e~ — > qq events than 7717^. The final event count 
is performed in the signal region 1.3 < m < 1.8 GeV/c 2 . 
Figure n illustrates the pseudo mass spectra of simulated 
signal and background contributions after the topology 
selection. 

For this analysis, a comparison of Monte Carlo simu- 
lation and data has shown that e + e~ — > qq background 
contributions cannot reliably be extracted from simula- 
tion due to difficulties in modeling the fragmentation 
processes. The shape of the pseudo mass distribution 
appears to be correctly modeled, but the overall normal- 
ization is not. Therefore, we estimate this background 
directly from data, using a method described below. On 
the other hand, the simulation of t pair events yields a 
reliable estimate of their expected background contribu- 
tion, verified by loosening requirements that suppress the 
tau background. A breakdown of the signal efficiencies 
and individual r background contributions are listed in 
Table|l]for each selection step. In this table the generic r 
sample does not contain 5ir and 57T7r° decays, which are 
listed separately since they comprise the only background 
from r decays after the final selection. 

The efficiencies for events in the signal region for t~ — > 
Att~3tt + is t and r~ -> 47r-37r+7rV are (9.4±0.1±0.6) % 
and (9.3 ± 0.1 ± 0.6) % (see Table 0. The first errors are 
statistical and the second systematic. 

The systematic errors on the signal efficiencies in- 
clude contributions from uncertainties in the reconstruc- 
tion of charged tracks (5.2%), the uncertainty associ- 
ated with the particle identification on the signal and 
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FIG. 1: Pseudo mass distributions of the seven charged tracks 
after the 1-7 topology selection for signal (dashed-dotted: 
r~ — » 47r~37r + 7r°i/ T ; dotted: r~ — > 47r - 37r + z/ T ), generic r 
(including 5-prong decays) (dashed), and e + e~ — > qq Monte 
Carlo events (solid). Entries are normalized to the data lumi- 
nosity. To display the two signal modes we assumed a branch- 
ing fraction of 2 x 10 -5 . 



tag side (2.7%), the luminosity measurement and the r 
pair cross-section determination (2.3%), and the uncer- 
tainty on the generic 1-prong r decay branching ratios 
(0.5%). Since the efficiencies of the two decay channels 
to the 7-prong r decays are in good agreement, we aver- 
age them and obtain an overall efficiency for the decay 
t~ -> 47r-37r+(7r K of (9.4 ± 0.6) % which includes the 
statistical and systematic contributions described above. 

To determine the number of e + e~ — > qq events in the 
signal region, we use the following procedure: we his- 
togram the pseudo mass distribution of all data events 
that satisfy the topology requirement and are in the 
pseudo mass sideband 1.8 < m < 2.6 GeV/c 2 . The con- 
tribution of the r background, determined from the simu- 
lation, is subtracted from this histogram, and the result- 
ing distribution is fit with a Gaussian function. We then 
fix the mean and the width of this Gaussian, and use it 
to fit the r-background-subtracted distributions of side- 
band events for all subsequent cuts, floating only the nor- 
malizations. Integrating the area of the resulting Gaus- 
sian functions in the signal region 1.3 < m < 1.8 GeV/c 2 
yields an estimate of the e + e~ — > qq background contri- 
butions. 

Tablelnlhsts the number of expected t and e + e~ — * qq 
background events in the signal region after the four se- 
lection steps. In Figure |21 (a) - (d) we show the tau- 
background-subtracted pseudo mass data distribution af- 
ter the four selection steps. Although entries below 
1.8 GeV/c 2 are shown here, these events were hidden dur- 
ing the development of our analysis to avoid experimenter 
bias. Overlaid in the individual figures are the fit re- 
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TABLE II: Predicted and observed number of events in the 
signal region of 1.3 < m < 1.8GeV/c 2 . The r background 
yield is obtained from the simulation, and its error reflects 
finite Monte Carlo statistics. The e + e~ —* qq yield is 
determined by fitting the data in the pseudo mass sideband, 
and its error results from data statistics. 





T Bg. 


e + e -> qq Bg. 


Observed 


1-7 Topology 


128 ± 13 


574 ± 21 


695 


Particle ID 


28 ±6 


241 ± 10 


244 


Conv. Veto 


2.4 ± 1.3 


119 ±5 


104 


1-Prong Tag 


1.3 ± 1.0 


20.3 ± 0.7 


24 



function, which convolves a Poisson distribution with two 
Gaussian resolution functions, accounting for the uncer- 
tainties in the background and in the efficiency: 



C(B;n,b,f,b,f) 



1 



TV. 27T(7b<7/ 

(2) 

where B denotes the branching fraction of 
t~ — > 47r _ 37r + (7r°)^ T , / = 2N TT e, e is the signal 
efficiency, /i = (n) = fB ± 6, with n the number of ob- 
served events, and b and at, are the expected background 
and its error, tjf incorporates the errors on the signal 
efficiency and the number of r pair events. 



TABLE III: Observed and predicted number of e + e~ — » qq 
background events in the 1-8 and 1-7 topology data for the 
different selection criteria. 





1-8 (fixed) 


1-8 (free) 


1-8 obs. 


1-7 (free) 


Topology 


19.0 ± 2.7 


19.0 ±2.7 


23 


574 ± 21 


Particle Id 


12.2 ± 1.6 


11.2 ±2.0 


10 


222 ± 19 


Conv. Veto 


2.7 ±0.3 


2.6 ± 1.4 


1 


126 ± 18 


1-Pr. Tag 


0.5 ±0.1 


0.3 ±0.5 





20.2 ±7.7 



suits to the pseudo mass spectra in the range of 1.8 to 
2.6 GeV/c 2 . 

After the final selection we observe 24 events in the 
data (see Figure |21 (a) ) with a total number of predicted 
background from r and e + e~ — > qq events of 21.6 ± 1.2. 
The statistical error of the e + e~ — > qq background esti- 
mate is derived from the statistical uncertainties of the 
parameters in the Gaussian fit. 

To validate the e + e _ — > qq background estimation we 
use 1-8 topology data 01 > which has negligible expected 
signal contributions. We fit the pseudo mass distribu- 
tion between 1.8 and 2.6 GeV/c 2 , integrate the fit func- 
tion in the pseudo mass region 1.3 to 1.8 GeV/c 2 , and 
compare this with the number of events found in data 
(see Table ITTTl columns two and four). In addition, we 
repeat the fit with mean and sigma floating ("free" fit) 
in each individual distribution to compare the expected 
numbers of events (Table ITTTl column three). As a fi- 
nal cross check, fits with the Gaussian mean and width 
floating were performed for the 1-7 topology data, and 
the results are included in Table IIIII column five. We 
note the very good agreement of the fits with fixed or 
floating Gaussian parameters. 

The systematic error on the number of expected back- 
ground events is ±0.4 events, derived from variations 
of the fit range and shape of the extrapolation func- 
tion. The total number of expected background events is 
21.6 ± 1.3. 

We calculate the branching fraction of the r~ — > 
47r _ 37r + (7r°)^ r decay based on the following likelihood 



The likelihood function is maximized with respect to 
the branching fraction B, f and b, and the following nu- 
merical value for the branching fraction is obtained by 
MINUIT0: 

B(j- -► 4^-3^+(7r°K) = (0.7±i;|) x 1(T 7 . 

Since we have no evidence for a signal we compute a 
Bayesian upper limit using a uniform prior in the branch- 
ing fraction, the background, and the efficiency. This is 
done by integrating out / and b in the likelihood function 
and plotting C as function of B. In this way we normal- 
ize the distribution to unity and get as the result of this 
analysis an upper limit at the point where the integral 
reaches 0.9: 

B(t~ -> 47r-37r + ( 7 r ( y T ) < 3.0 x 10~ 7 (at 90% CL). 

With the same approach, setting the number 
of observed events N f, s to the expected num- 
ber of background events of N exp = 21.6, we 
calculate the sensitivity of the analysis to be 
B N "»= N ""(r- -> 4ir~3ir+(ir )v T ) < 2.5 x 10~ 7 at 
90 % CL. 

In addition to this inclusive result, we set limits on 
the branching fractions of the exclusive decay modes 
r~ — > 47T~37r + iA r and r~ — » 47r~37r + 7r°^ T . To select 
r~ — » 47r~37r + ^ 7 - candidates in the inclusive sample, we 
require the number of photons on the signal side to be 
zero. This yields a signal efficiency for r~ — > 47r~37r + z/ r 
decays of (5.5 ± 0.3) % while reducing the expected 
generic r decay and e+e~ — > qq background to 3.9 ± 0.8 
events. The decays t~ — > 47r - 37r + 7r ( V r are treated 
as background in this case and have a reconstruction 
efficiency of (0.8 ± 0.1) %. Reversing this selection by 
demanding at least one reconstructed ir° on the signal 
side yields a signal efficiency for the r~ — » 47r~37r + 7r°z; T 
decay of (3.6 ± 0.3) % and an expected generic r decay 
and e + e~ — > qq background of 8.2 ± 0.5 events. In this 
case, the reconstruction efficiency for t~ — > 47r~37r + i/ T 
is (0.3 ± 0.0) %. The systematic uncertainties of the 
r~ — > 47r _ 37r + ^ T mode are identical to the inclu- 
sive measurement already discussed above. For the 
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Pseudo Mass (GeV/c 2 ) 

FIG. 2: Pseudo mass distributions for 1-7 topology data, (a) 
after conversion veto; (d) after 1-prong tag. 

t~ — > 47r - 37r + 7r°M r mode, an additional uncertainty of 
5.0% on the efficiency of the ir° reconstruction is taken 
into account. 

The likelihood function (|5J) was modified to accommo- 
date one exclusive mode acting as background for the 
other. We observe eight events in the r~ — > 4:Tt~3tt + v t 
signal region and seven events in r~ — > 4ir~ 3ir + ir a v T . We 
calculate the following numerical values for the branching 
fractions: 

B(t~ -v 47r"37r+^ T ) = (2.0±i;|) x 10~ 7 
S(t" -» 47r"3vr + 7r ^ r ) = (-1.0 ± 1.8) x KT 7 . 

Without evidence for a signal we compute Bayesian 
upper limits using uniform priors in the branching frac- 
tions, the backgrounds, and the efficiencies, of 

B(t~ -> <lTr-3Tr + v T ) < 4.3 x 10~ 7 (at 90% CL) 
B(r~ -> 47r"37r+7rV T ) < 2.5 x 10~ 7 (at 90% CL). 

With the same approach, setting the number of ob- 
served events N a b s to the expected number of background 
events, we calculate the sensitivities B oba ex p(r~ — > 
47r-37r+i/ T ) < 2.2 x 10~ 7 and B N -""= N ^" (t~ -> 
47r-37r+vr ^ T ) < 4.2 x 10~ 7 . Figures (b) and (c) show 
details of the data pseudo mass spectra with an overlay 
of the expected background distributions. 




after 1-7 topology selection; (b) after particle identification; (c) 



This analysis improves the existing experimental lim- 
its by an order of magnitude for the inclusive mode, but 
is still several orders of magnitude larger than the theo- 
retical prediction. The exclusive decays are reported for 
the first time and are consistent with the inclusive result. 
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FIG. 3: Details of the data pseudo mass spectra, (a) r 
candidates. 



4ty 37r + (7r°)^T, (b) r — > 47r 3tt + v t , (c) r — > 47r 37r + 7r l V T event 
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